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We examine the dynamics of the electroweak phase transition in the early 
Universe. For Higgs masses in the range 46 < Mh < 150 GeV and top quark 
masses less than 200 GeV, regions of symmetric and asymmetric vacuum co- 
exist to below the critical temperature, with thermal equilibrium between the 
two phases maintained by fluctuations of both phases. We propose that the 
transition to the asymmetric vacuum is completed by percolation of these 
sub-critical fluctuations. Our results are relevant to scenarios of baryogenesis 
that invoke a weakly first-order phase transition at the electroweak scale. 


PACS numbers: 98.80-k, 98.80.Cq, 12.15.-y 
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The realization that gauge symmetries can be restored at high temperatures, com- 
bined with the success of the big-bang model of cosmology, has generated a lot of interest 
in the study of cosmological phase transitions. 1 First-order phase transitions are char- 
acterized by an energy barrier separating the symmetric from the asymmetric phase at 
the critical temperature Tc when the two phases have equal free energy. First-order 
transitions may generate out-of-equilibrium conditions, which can have important effects 
upon the properties and evolution of the early Universe. Two well-known examples are 
models of inflationary cosmology that invoke a first-order transition at the grand-unified 
scale, 2 and the production of inhomogeneities at the quark-hadron transition. 3 

In this letter we study the electroweak phase transition in the minimal (i.e., one Higgs 
doublet) model. We will restrict our study to fairly light Higgs masses, ranging from 46 
GeV up to 150 GeV, and top-quark masses in the range 100 to 200 GeV. For this range 
of parameters the phase transition is weakly first-order. The high temperature minimum 
of the potential is the symmetric state (<f>) = 0. At some temperature T\ > Tc the 
potential develops a local asymmetric minimum at <f>+ > 0. As the system cools, the 
difference in free energy between the symmetric and asymmetric state decreases; finally 
at the critical temperature Tc the asymmetric minimum is degenerate with the symmetric 
minimum. Below Tc the asymmetric minimum has the lower free energy. Eventually, at 
some temperature T 2 < Tc the symmetric minimum becomes unstable. 

Two scenarios have been proposed for the completion of such transitions. In the 
“standard” picture, the Universe remains in a homogeneous state of symmetric vacuum 
below Tc, until the symmetric state becomes unstable at Tj. Then the field evolves 
classically to the asymmetric minimum. 1 Recently a second scenario has been proposed 
where again the Universe remains in a homogeneous state of symmetric minimum to Tc, 
then between Tc and T 2 the homogeneous state is terminated by nucleation of bubbles 
of asymmetric (true- vacuum) phase which grow and eventually percolate the volume. 4 
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We propose that the transition is completed by a new mechanism: percolation of sub- 
critical fluctuations of the asymmetric phase. We argue that by the time the Universe has 
cooled to Tc, the vacuum is not a homogeneous state of symmetric vacuum, but rather 
an emulsion of symmetric and asymmetric vacua, each existing with equal probability. 
Below Tc the fraction of the Universe in the symmetric state gradually decreases, and 
the transition is completed by percolation of many regions of asymmetric phase. 

We use a method developed by Gleiser, Kolb, and Watkins (GKW) designed to study 
the approach and maintenance of thermal equilibrium in phase transitions. 5 In this 
approach the thermal fluctuations of the field are modeled by the creation of regions 
(bubbles) of one phase inside of the other. These fluctuation regions are spherical and 
have a size of the thermal correlation length of the Higgs field, L GKW use detailed 
balance to find the rate of creation of fluctuation regions of false vacuum inside a true- 
vacuum region to be i~ x exp(— Af’/T) where A F is the difference in free energy of the 
region and the homogeneous state. If these rates are large compared to the expansion 
rate H , then the relative population of the phases should be distributed according to 
Boltzmann statistics. We find this to be the case for the electroweak transition with 
top and Higgs masses in the aforementioned ranges. Our results should be relevant to 
the recently proposed scenarios of baryogenesis at the electroweak scale, which naturally 
invoke out-of-equilibrium conditions during a first-order phase transition. 6 

In the study of phase transitions the Higgs field (or its equivalent) plays the role of 
the order parameter. In practice, when the system is initially in thermal equilibrium, 
the study of the phase transition reduces to the construction of the finite-temperature 
1-loop effective potential, which incorporates the interactions of the Higgs field with itself 
and with other fields in the model at some temperature T. 7 The effective potential is 
equivalent to the homogeneous part of the free energy and its minima determine the 
equilibrium properties of the system. We neglect contributions of the Higgs field to the 
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1-loop potential. This should be valid for Higgs masses below about 150 GeV. In order 
for the potential to be stable with these small Higgs masses, the top quark must be less 
than about 200 GeV. 8 The transition can be studied using a high-temperature expansion 
of the effective potential which, as shown by Turok and Zadrozny 9 and by Anderson and 
Hall, 4 is very reliable in the relevant range of temperatures. They obtain for the potential 
(we will follow the notation of Ref. 4) 

V(<f>,T) = D (t 2 — T 2 ) <f> 2 - ET<f> 3 + (1) 

where the constants D and E are given by D = [6(mw/<r) 2 + 3 (m^/cr) 2 + 6 (mx/<j) 2 ] /24, 
and E = [6 (tow/g’) 3 + 3 {mzl<r) 3 \ /127r. Here T 2 is the temperature at which the origin 
becomes an inflection point (i.e., below T 2 the symmetric phase is unstable and the field 
can classically evolve to the asymmetric phase), and is given by Ti = — 8f?<r 2 )/4Z?, 

where the physical Higgs mass is given in terms of the 1-loop corrected A as m 2 H = 
(2A + 12B) <r 2 , with B = (6m^- + 3m| — 12my) /647r 2 <r 4 . We use mw = 80.6 GeV, 
mz = 91.2 GeV, and <r = 246 GeV. The temperature-corrected Higgs self-coupling is 

‘in (mJ./c f T ! ) , (2) 

where the sum is performed over bosons and fermions (in our case only the top quark) 
with their respective degrees of freedom 9b{F)i and In eg = 5.41 and lncp = 2.64. 

Apart from T 2 , there will be two temperatures of interest in the study of the phase 
transition. For high temperatures, the system will be in the symmetric phase with 
the potential only exhibiting one minimum at ( <f> } = 0. As the Universe expands and 
cools an inflection point will develop away from the origin at <f> = 3 ET x /2Xt, where 
T\ is given by Ti = Ti/yjl — 9E^/BXtD. For T < T 1? the inflection point sepa- 
rates into a local maximum at <j>- and a local minimum at <f >+ , with <f>± = {3 ET ± 
[I *E 3 T 2 - 8A t D{T 2 - T 2 )] 1/3 }/2A r . At the critical temperature T c = T 2 f y/l -E 2 /\ T D, 
the minima have the same free energy, V(<J> + ) = V(0). 


Ar = A - i-L K/csT’) + X> (2E) 
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In the usual picture of a first-order transition, the field starts in thermal equilibrium 
in its symmetric minimum at (<f>) = 0, and as the Universe cools below Tc the symmetric 
phase becomes metastable and decays by nucleation of bubbles of the asymmetric phase: 
bubbles of size greater than the critical size grow, converting the symmetric phase into 
the asymmetric phase. The success of this scenario depends crucially on the assumption 
that the field is in a homogeneous state of the symmetric minimum as the Universe 
cools below Tc- However, hot systems tend to fluctuate around their equilibrium states, 
and the probability to find the system in a state other than its ground state has a 
relative probability given by the Boltzmann factor, exp [— F(T)/T\, where F(T ) is the 
free energy for the particular fluctuation. For high enough temperatures and slow enough 
cooling rates, the system will have a large probability to populate other accessible states. 
For a system with a metastable and a true-vacuum state the equilibrium probability is 
exp [— (A F(T)) /T], with A F(T) being the free energy difference between the two states. 
For the electroweak model with the potential given by Eq. (1), as the temperature drops 
below T\ thermal fluctuations may drive the system into equilibrium populating the new 
minimum at <f> + . If this is the case, as the temperature drops below Tc the Universe will 
be filled by a two-phase emulsion, and the kinetics of the transition will be quite different 
than the usual false vacuum decay scenario. 

GKW assume that the dominant statistical fluctuations are sub-critical bubbles of 
roughly a correlation volume which interpolate between the two minima of the free en- 
ergy. Denoting the minima for the electroweak model <f> 0 and 4> + , for the symmetric and 
asymmetric states, the rates for fluctuations between the two states are 

rfrjjo-t+j 2 ; mo(T)exp [— ,F + (T)/T] ; r(3') (+ _ <1 , ~ m„(T)exp hF 0 (T)/T ) , (3) 

for a fluctuation of the asymmetric (symmetric) phase within a region of the symmetric 
(asymmetric) phase. F+(T) is the free energy of a fluctuation of the asymmetric phase 
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and Fo(T ) is the free energy of a fluctuation of the symmetric phase. For simplicity, 


we assumed the same correlation length [^(T) -1 = to 0 (T) = y/V"(0, T) j around the two 
minima. Now we must estimate the free-energies F + (T ) and F 0 (T). The free energy of a 
fluctuation in the order parameter is given by (for details see GKW) 

F(T) = J cPx [i { V^) ! + V(*. r)] , (4) 

where V{<j> , T) is given by Eq. (1) and the order parameter tf> is the amplitude of the Higgs 
field. We are interested in fluctuations of roughly a correlation volume that convert re- 
gions of symmetric phase into regions of asymmetric phase and vice-versa, which will give 
the dominant contribution to the transition amplitude. Since these field configurations 
are not solutions of the euclidean equations of motion, we adopt a variational approach 
to determine the dominant configurations with minimal free energy. Thus, we take for 
the sub-critical bubbles, 

<f>+{r) = <£+exp (“ r V* 2 ) i M r ) = 0+ I 1 ~ exp (-r 2 /^ 2 )] , (5) 


where 0 + ( O )(r) is an 0(3)-symmetric bubble of asymmetric (symmetric) phase nucleated 
in the symmetric (asymmetric) phase. Introducing the dimensionless variables X(p) = 
<f>(r)/<r , £(T) = ^(T)<r, 0 = T/<t, and p = r<r, the free-energies are given by 


f+W = [M + p (^ («’ - «!) - 


and 


Fo(0) = T 3/2 Xll<r 


TT + p (f (<’ - »l) (V5 - 8) 


( 6 ) 


(-1 + ^ + i) + E 6 (3 - ^ x + ) . ( 7 ) 

The free energy F+(T)/T for T — Tc is shown in Fig. 1 as a function of the Higgs 
mass for several values of the top mass. This free energy will determine the equilibration 
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properties of the system as the temperature drops below T\. Note that F+(T) increases 
as the temperature drops. This is a consequence of the fact that the free energy is 
dominated by the gradient energy, and as the temperature decreases the asymmetric 
minimum moves away from the origin. In order to establish thermal equilibrium by 
overcoming the energy barrier, the thermal fluctuation rate in going from <j> = 0 to 
<f> = <f> + must be large compared to the expansion rate of the Universe: r[o-.+]/-£f 1 } 

with H ~ 1.66 gV 2 T 2 /Mpi, and g . ~ 110 is the number of effective relativistic degrees of 
freedom at the electroweak scale. Negleting pre-factors, this condition can be easily seen 
to lead to the inequality F + (T)/T 34. 

From Fig. 1 we see that for most of the parameter space studied F+(Tc)/Tc is com- 
fortably less than the critical value of 34, so equilibrium should be established at Tc by 
fluctuations going from (<f>) = 0 to ( <f> ) = 0+. 10 (Obviously, fluctuations in the opposite 
direction will have smaller free energy until T = Tc when the two free-energies are the 
same. 5 ) Thus, we conclude that at Tc the Universe is not in a homogeneous state of 
symmetric vacuum as assumed in all previous works on the subject. Another indication 
that large fluctuations in . the Higgs field will be important is to note that at T = Tc, 
<t>+/T = 2E/Xt ~ 2 x 10“ 2 /A t* Of course At depends upon mj and T, but using its 
tree-level value of Ao = nig /2a 2 = 0.04(rrifl-/100 GeV) 2 is a good approximation 4 and 
shows that <f> + /T at Tc is never much greater than unity, and typically is less than unity. 
Since T sets the scale for thermal fluctuations, the system should “feel” both minima. 11 

So far we have established that (for Higgs and top-quark masses we considered), it is 
quite easy for equilibrium of the two vacuum states to be achieved. In this case as T drops 
below Tc, the Universe will be filled by a two-phase emulsion, with rapidly fluctuating 
regions of symmetric and asymmetric phases, separated by roughly a correlation volume. 
As T drops below Tc fluctuations from the asymmetric phase back to the symmetric phase 
become more and more suppressed, and the asymmetric phase will occupy more than 50% 
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of the Universe. The mechanism by which the transition is completed is complicated and 
will depend on the temperature at which the fluctuation rate freezes out, TV- If TV > T 2 , 
the symmetric phase is still locally stable, and correlation volume regions of this phase 
will shrink under surface tension, while regions of the asymmetric phase, having lower 
free energy will percolate. In Fig. 2 we show the ratio of both rates to the expansion 
rate as a function of the temperature for mg = 60 GeV and mx = 130 GeV. Only for 
fairly light Higgs will 2> be larger than T 2 . For TV < T 2 , the symmetric phase becomes 
unstable and fluctuations to the symmetric phase can classically roll back down to the 
asymmetric phase. The Universe will be quickly permeated by the asymmetric phase, 
since any interface region is energetically disfavored and will move toward the symmetric 
phase converting it into the true vacuum. 

We have shown that for the minimal standard model, with 46 < m# < 150 GeV and 
my < 200 GeV thermal equilibrium will be maintained during the electroweak phase 
transition despite the fact that there is an energy barrier between the two phases. We 
assumed that equilibrium is maintained by the thermal nucleation of sub-critical field 
configurations of roughly a correlation volume, since these are the statistically dominant 
fluctuations at temperature T. The free energy of these configurations was estimated by 
assuming they are 0(Z) symmetric and that they interpolate between the two phases. 
Preliminary numerical simulations indicate that the ansatz used here is correct within 
10%, although a more detailed analysis is necessary. In any case, the lesson is clear; 
given enough time and heat, a system will thermalize due to the nucleation of field 
configurations that convert one vacuum into another. This result can be easily extended 
to the recently proposed scenarios of baryogenesis at the electroweak scale. Of course, 
we must go beyond the standard model since it does not have a source of CP violation, 
and extensions keeping only a Higgs doublet 12 and with two Higgs doublets 13 have been 
proposed. A successful baryogenesis scenario cannot assume a metastable symmetric 
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phase below Tc only because there is a barrier between the two phases. 

We would like to thank Scott Dodelson, Erick Weinberg, and Michael Turner for very 
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FIG. 1. The free energy of the sub-critical fluctuation at the critical temperature as 
a function of the Higgs mass for several values of the top-quark mass. 
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FIG. 2. The ratio of the fluctuation rate to the expansion rate as a function of tem- 
perature. 
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